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Abstract. We study the evolution of the orbit of a planetary companion as the primary star evolves 
during the red giant branch phase. Our aim is to determine a prescription for the capture radius and 
capture timescales for companions in a wide range of masses, from planet to star, including the 
effect of a non constant / parameter. Our initial results reproduce past findings. According to our 
study the mass-loss during the red giant branch does not significantly affect the semi major axis of 
the planet and for a 1.5 M star the capture radius is approximately 300 R or 1.4 AU. 
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INTRODUCTION 

The evolution of a star with a planetary companion does not differ a lot from the 
evolution of a single star until the planetary companion becomes too close to the star 
due to stellar expansion or orbital evolution. If the companion gets too close to the star it 
can be tidally captured. The farther out a companion can be captured by an evolving star, 
the larger the birthrate of stars that derive from binary and planetary interactions (e.g., 
cataclysmic variable stars). Population synthesis models [6] need to treat tidal capture if 
they are to make credible predictions of binary population densities and birthrates (e.g., 
the type la supernovae). 

Though many studies have been conducted on the subject (e.g., Villaver and Livio 
9), very few have taken into account the effect of synchronisation on the evolution of 
the orbital parameters. In the case of planetary companions, synchronisation does not 
affect the orbital evolution, but when companions get heavier, as in the case of brown 
dwarfs and low mass main sequence stars, synchronisation can play a major role on the 
evolution of the orbit. When an object starts to tidally synchronise, it slows the rate of 
the evolution of the orbital separation. In the absence of mass loss, the orbit will not 
evolve in a perfectly synchronised system. The long term purpose of this study is to 
find a general prescription for tides on orbits of planets and stellar companions that goes 
beyond what has been done by [7]. 

In addition, our study will encompass the effect of the dimensionless parameter / [5], 
which depends on the convective mixing length parameter. For a larger value of / the 
strength of the tides will be stronger and the evolution of the orbit will be faster. Most 



studies use / = 1 (e.g., Verbunt and Phinney 8, Soker 7), but [3] has shown that / need 
not be constant to match observations of post-main sequence binaries, which means that 
the functional dependence of / remains unconstrained. 



The evolution of the semi major axis is governed by competing processes. During the 
red giant branch (RGB) phase of the star, tidal force, gravitational and frictional drag, 
variation in mass of the star, and changes in the planet mass due to accretion of stellar 
gas as well as ablation by the star, all play their parts. However, e.g., [9] established 
that during the RGB phase, gravitational and frictional drag, along with the mass gained 
by the planet by accretion do not contribute significantly to the evolution of the semi 
major axis. Furthermore, in accordance with the study of [3], we assume the mass of 
the secondary remains constant through the evolution. For this study we consider alMj 
planet revolving around a 1.5 M star, with zero eccentricity. 

The orbital change can be expressed as the sum of the change due to mass-loss, and 
that due to tidal interaction 



where a is the semi major axis, M* is the mass of the star and d and M* are their time 
derivatives (e.g., Nordhaus et al. 3). 

The tidal force plays a very important role in the evolution of the semi major axis 
during the RGB phase of the star. For RGB stars, the surface rotation period is much 
less than the orbital period of the secondary. In this situation, the sign of the tidal torque 
is such that the orbit decays. Even though there are various tidal prescriptions discussed 
(e.g., Hut 2), the most widely used for convective stars is that formulated by [12, 13, 11]. 
Tidal decay in this mechanism is governed by convective eddies in the star's envelope: 



where (0 is orbital frequency of the companion, £1 is spin frequency of the star, M env is 
the mass of the envelope, R* is the radius of the star and q is the mass ratio between 
the planet and the star. We are also including the effect of synchronisation between 
the orbital frequency of the companion and the spin frequency of the star, though it 
plays very little role in the evolution of the orbit for a planetary mass companion. The 

parameter / is a dimensionless quantity, equal to 1 .01 ( ") 3 , where a is the mixing length 
parameter [8]. For the purpose of the current study we have maintained / constant, 
although it will be varied in the future. The parameter x is the convective eddy turnover 
timescale in the envelope. 

We are using the stellar evolution code, MESA of [4] for the evolution of the star. We 
evolved a binary system consisting of a 1.5 M star and a companion of 1 Mj through 
the entire RGB phase of the star. During each time step equaling 100 years, we update 



EVOLUTION OF THE ORBITAL SEPARATION 




(1) 




the value of the stellar parameters and calculate the new value of the a in accordance 
with [11]. The parameter a is taken to be 2 in accordance with [10].The initial spin of 
the star at the begin of the RGB is taken to be 3 km s 1 [1]. 

RESULTS 

We arbitrarily start our code when the primary has a radius of 5 R . We run three tests 
where we vary the initial separation of the 1 Mj : 275 R , 320 R , 360 R . Our 
preliminary findings are similar to those of [9] and [3]. We observe that the mass- 
loss during the RGB does not contribute as much to the evolution of the semi major 
axis, a, of the binary even though the 1 .5 M star has lost roughly 8 percent of its initial 
mass. Figure 1 shows the evolution for the semi major axis for entire RGB phase of 
the star. The increase in a towards the end of the RGB evolution is due to mass-loss. 
The mass-loss during the AGB will have a much larger effect. In addition, the detailed 
calculations of the MESA code during the thermally pulsating AGB are likely to result 
in a different evolution of the orbital separation compared to tidal codes that use off the 
shelf stellar evolutionary tracks (e.g., Nordhaus et al. 3). Figure 2 just focuses on the 
last 200,000 years of RGB evolution. The two plots show a range of initial separations 
for which the companion may get captured or evade capture. The thick continuous line 
shows the evolution of the primary's radius, the three discontinuous lines describe the 
evolution of the semi major axis a. As the system evolves we find that the distance 
beyond which the companion does not get captured is approximately 320 R or 1.5 AU. 
This is approximately 2.5 times the maximum RGB radius of this star. Further study 
needs to be done in the cases where the companion is in the range 300 - 320 R , here 
the companions may get captured during the beginning of the horizontal branch phase 
of the star. 



Future Work 

We expect to find that since mass-loss dominates during the AGB phase of the star, 
it is unlikely that the companion will get engulfed if it is much farther than a couple 
of stellar radii. This was also the result of [3]. We will extend the work of [3] to more 
massive companions, including the effect of synchronisation. By including a variable 
/ value we will extend the results of [9], and by including a proper treatment of the 
thermally pulsating AGB we will have a better representation of the AGB star radial 
evolution, determining when and how companion get engulfed during this phase. 
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FIGURE 1. Evolution of the radius of the star during the entire RGB phase (continuous line) and the 
orbital evolution of the semi major axis with three initial separations (275, 320, 360 R©, discontinuous 
lines). 
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FIGURE 2. Evolution of the radius of the star during the last 200,000 years of the RGB (continuous 
line) and of the orbital evolution of the semi major axis for 3 distinct initial separations during the same 
period (275, 320, 360 R , discontinuous lines). 
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